Relativistic Scalar Gravity: 
A Laboratory for Numerical Relativityo 
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We present here a relativistic theory of gravity in which the spacetime metric is derived from a 
single scalar field $. The field equation, derived from a simple variational principle, is a non-linear 
flat-space four-dimensional wave equation which is particularly suited for numerical evolution. We 
demonstrate that while this theory does not generate results which are exactly identical quantita- 
tively to those of general relativity (GR), many of the qualitative features of the full GR theory are 
reproduced to a reasonable approximation. The advantage of this formulation lies in the fact that 
3D numerical grids can be numerically evolved in minutes or hours instead of the days and weeks 
required by GR, thus drastically reducing the development time of new relativistic hydrodynamical 
codes. Scalar gravity therefore serves as a meaningful testbed for the development of larger routines 
destined for use under the full theory of general relativity. 



PACS numbers: 04.25.Dm, 04.30.Db,04.40.-b 



I. SCALAR THEORIES OF GRAVITY 

Computational loads in numerical relativity are typi- 
cally immense. Fox [nl estimates that the computation of 
the gravitational waves from black hole collisions would 
require approximately 100,000 hours on a Cray Y-MP. 
Current [Q] improvements in parallel machines and algo- 
rithms may cut this to 50 hours on rare machines. It is 
very difficult to do even test problems under these con- 
ditions. Clearly, one cannot wait days or weeks to get 
results to use in deciphering the cause (or even detecting 
the existence) of some problem in the code. When one 
considers that we are not even completely certain what 
new problems we will be forced to face, the situation is 
completely unacceptable. While many problems will be 
unique to the Einstein equations and must be faced in 
that context, others, notably in relativistic hydrodynam- 
ics, will be similar using any metric theory of gravity. 

Shapiro and Teukolsky |^ have proposed that instead 
of starting with the full equations of general relativity, 
we should instead look at a simpler relativistic theory of 
gravity in order to discover and learn how to resolve these 
unknown issues while utilizing far less computational re- 
sources than the full problem would require. The scalar 
theory presented here is intended to provide a testbed 
for relativistic hydrodynamics. It is to be implemented 
in code that accepts a matter stress-energy tensor and 
generates an updated metric at each time step. The hy- 
drodynamic code that accepts this metric and updates 
the hydrodynamic variables and the stress-energy ten- 
sor can be exactly the same code that is to interact in 
the same way with a program that solves the Einstein 
equations. Thus speeded development and testing of hy- 
drodynamical codes is the principal motivation for de- 
veloping this scalar theory. Any gravitational physics or 



techniques that may be carried over to Einstein equation 
solvers will be a bonus. 

Shapiro and Teukolsky's prescription is to explore the 
scalar gravitation theory presented in Exercise 7.1 of Mis- 
ner, Thorne, and Wheeler (MTW) @]. This theory, how- 
ever, violently disagrees with observations and fails all 
three of the classical tests of general relativity. It is not 
clear that this theory is close enough to general relativity 
to prove useful to our purposes. Nevertheless, the poten- 
tial benefits of a simplified testbed are such that a further 
examination of scalar gravity theories is warranted. 

Ni ||5[] provides a compendium of metric theories orga- 
nized in terms of their "Parameterized Post-Newtonian 
limits" as defined by Thorne and Will Q. In general, 
scalar theories can be grouped into three classifications: 

1. Conformally flat theories 

2. Stratified conformally flat (SCF) theories 

3. Other scalar theories 

Conformally flat theories are some sense the simplest pos- 
sible theories of gravity. In this class of theories, the 
metric takes the form 



rfs^ = ^{-df + dx^ + dy^ + dz^) 



(1.1) 



where the variable ^ is some function of the scalar field 
<&, e.g., '^ = exp2$. We use the symbol $ through- 
out for the basic scalar gravitational field, normalized so 
that in the Newtonian limit it becomes the Newtonian 
gravitational potential, asymptotically $ = ~M/r = 
—GM/rc^. We shall also use units with G = 1 = c where 
G is the Newton-Cavendish gravitational constant. The 
field equation (which determines <&) can technically be 
anything, but is generally some variation of the scalar 



wave equation, with some version of the matter density 
serving as the source term. Ni describes five different 
examples of this type of theory. [|| 

All conformally flat theories (including the theory of 
MTW Ex.7.f) suffer from a fatal flaw: since MaxweU's 
equations are conformally invariant, all conformally flat 
theories must predict zero light deflection, in obvious 
disagreement with experiment. Most also fail various 
other solar-system tests of general relativity, but the zero 
light-deflection prediction is a fundamental consequence 
of each of these theories and cannot be avoided. 

One way to solve the light-deflection problem is to pos- 
tulate that there exists a universal preferred reference 
frame and that while spacelike slices of this spacetime 
(referred to as "strata") are conformally flat, the full 
spacetime is not. The general form for the metric of 
this "stratified conformally flat theory" is 



ds^ = -^idf + *2(rfa;^ + dy^ + dz^ 



(f.2) 



where again ^i and ^2 are functions of the scalar field 
$. This type of metric obviously destroys Lorentz in- 
variance, however this is not as strict a requirement as 
one might think. Global Lorentz invariance is a treasured 
philosophical feature of general relativity with asymptot- 
ically flat boundary conditions, but in reality we usually 
model sources as though they existed alone and in isola- 
tion in the Universe. The fact that we have a preferred 
reference frame in this case is not such an important is- 
sue. 

The last category, more general scalar theories, has 
seen little examination in the literature and we do not 
study it here. 



II. SIMPLE RELATIVISTIC SCALAR GRAVITY 

A. Variational Principle and Metric 

We propose a variational principle which leads to a 
simplified theory closely related to Ni's "Lagrangian- 
based stratified theory" , which until 1972 was thought 
to agree with all tests of general relativity. It has since 
been found to disagree with the Earth rotation rate ex- 
periment |Tq] and with white dwarf stability observations 
M (although Ni notes that the validity of the latter ex- 
periment can be questioned since it is based upon ideas 
about pulsation-damping forces in white dwarfs which 
have yet t o be tested experimentally). As will be seen 
in Section HE , the PPN form of the simplified metric 
is identical to that found from Ni's theory, hence it sat- 
isfies (and fails) the same tests of general relativity as 
does Ni's theory. This stratified conformally-flat (SCF) 
scalar theory, therefore, has all of the advantages of Ni's 
more complicated theory, yet has a simpler form which 
eliminates many of the disadvantages. Additionally, the 
resulting field equation has the form of the non-linear 
Klein-Gordon equation, for which a significant amount 



literature exists, thus aiding in the numerical implemen- 
tation of the theory. 

We begin by selecting the Lagrangian: 



/ = -^ / ?/'"'*,,,*,.\/^rf^a; + / Anattcrrf^ 



(2.1) 



This differs from Ni's theory (cf. ||5|] or a summary in 
MTW ||, Box 39.1]) only in using r]^"' where Ni uses 
5'^" in the first term, so that our scalar field $ propa- 
gates along flat spacetime light cones. Here 77'^" is the 
Minkowski metric, and ^/—rj — 1 unless one is using 
nonstandard, e.g., spherical, coordinates. We generate 
the metric (with which all nongravitational fields inter- 
act just as in general relativity) in the same way as Ni: 



ds-" 



-e^^df+e-^^idx^ + dy^ + dz^) 



(2.2) 



We have briefly explored alternative forms of the strati- 
fied metric, 

ds^ = -(1 - -$)-"di2 _^M_ l<^\nr^^2 ^ ^y2 _^ ^^2n 

n n 

(2.3) 

but found that the most realistic (i.e., closest to Einstein) 
results are obtained when n is allowed to go to infinity. 



B. The Field Equation 

The field equation follows from varying $ in the action 
21) which gives 



dgfj,^ 






(2.4) 



The right hand side of this equation can be expressed in 
terms of the stress-energy tensor by using the definition 
found in MTW |, Chap. 21]): 



<s/:n 



^3lLV 



(2.5) 



The field equation derived from this variational principle, 
which is valid for a metric generated in any way from $, 
is therefore 



l~^-]U^ = -27r\/^ 



9$ 



(2.6) 



where y/—7]0 = d^^—'rjrj^^'^d,y. Utilizing the SCF metric 
~2h, this field equation becomes 



-f,U^^ 47r(rO° + e-'**%T*^') 



(2.7) 



The theory pro posed by Ni [pj differs only in using 
the metric (2.2) where our simplified theory uses the 
Minkowski metric 77^^ on the left hand side. 



C. Conserved Quantities 

It is often important to know which physical quanti- 
ties are conserved in a theory. These quantities are of 
particular importance if the field equation is to be solved 
numerically, as they are often used as to monitor the er- 
ror or stability of a numerical integration algorithm (e.g. 
Pi). A conserved energy and conserved linear and angu- 
lar momentum can be derived from the Lagrangian im- 
plicit in the variational principle of (2.1) using Noether's 
Theorem, since this Lagrangian density is invariant un- 
der spacetime translations and spatial rotations. Exposi- 
tions of Noether's theorems can be found, e.g., in Hill B] 
and in Byers ||l^ . For the present application the special 
case using the canonical stress-energy tensor suffices. For 
a Lagrangian density >C(^*, ^^,;i) which has no explicit 
coordinate dependence but may contain several fields ^ , 
the equations 



Only the last step in this equation uses the special SCF 
form ( |2.2| ) of g^i^ in terms of $. 

Since the i]'^'^ are constants, the conservation law 
Tx^ .y = can also be written T^'^ ^^ = by defining 
7"^7 = ij^'^Ta'^ . One then defines the total energy and 
momentum of the system as usual by 



pM 



^m0^3 



(2.15) 



and recognizes T'^'' as th e Poynting-like energy flux. The 
zero-component of (2.15) is the energy, defined as 



E = P° = - %°d 



Txl = 



where 



C5l 



dC 



-*^ 



,A 



(2.8) 



(2.9) 



with the corresponding energy conservation law 
^--fr da, . 



(2.16) 



(2.17) 



are easily shown to be a consequence of the field equations 



= 



a*A 



a. 



a*A, 



(2.10) 



Here T*"' = —% = S'^ is an energy flux like the 
Foynting vector. In most applications the bounding sur- 
face will be chosen outside the region occupied by mat- 
ter so only the gravitational terms remain as S^^.^^^ = 

-(cV4^G)$,o*,fc. 

Angular momentum conservation is a consequence of 
the symmetry of the stress-energy tensor. The gravita- 
tional part of T^'' is symmetric as a consequence of the 
Lorentz invariance of the Lagrangian >Cgiav from which it 



(See, e.g., pi|-|l3[.) Since the stress energy tensor Tx'' as is derived, and is also clear from its explicit deflnition 



defined in equation ( 2^ ) is linear in C = C 



-'grav 



A, 



it is also a sum of two parts and the energy-momentum 
conservation law can be written as 



1 

4Tr 



Cav) = J-V^l{-'^V''-'^,,<^..S''^ +V'^<i>,,^,.v'' 



[n 



Here 



"^(gr 



, 'y 



av) A 



(grav) A 
1 

47r^ 



-'(matter) A J -, 







^rj{~^^7r<P,,,<P,,6l 



(2.11) 



$,A*,p??^^) 



(2.18) 



However the matter Lagrangian is not Lorentz invariant 
in view of the preferred rest frame in which equation ( |2.2| ) 
holds. The matter stress-energy is defined by 



(2.12) 



-* (matter) 



^v^^gpc- 



sc 



— '/ yfSav y -t 



r/aj 



follows from equation (2.9) since £grav depends on only 
the single field $. Although £niattcr may depend on sev- 
eral matter fields, in it $ occurs only indirectly via g^^. 
Since we assume that iZmatter is generally covariant, one 
knows that its stress-energy tensor can equivalently be 
obtained as 



(2.19) 



and does not inherit the symmetry of T"^ except in its 
spatial components: T^mattcr)*'' = '?(mattcr)^' = e"^*T'^ . 
An angular momentum density can be defined as usual 



Mf""^ = x'^r^'^-x'^r^^ which satisfies Mf""^ . 



T'l^M. 



'?(matter)A — 2(7Aq 



sc„ 



Sga^ 



gxaV-gT 



(2.13) 



T^'^ so that a conserved angular momentum is defined 



To directly verify the con servation law (2.11) one uses the 
scalar field equation (2.7) and the covariant conservation 



J^3 



M'^° Sx 



(2.20) 



law for the matter fields T^'' -^ = in the form 



q- 7 FT ^^7 

-'(matter) A ^ -^ a A7 -* (matter) 



29a-i,\V-gT° 



= -$,A(r'" + e-*'^%rJ) 



(2.14) 



D. The Newtonian Limit 

In order to find the Newtonian limit for this theory, as- 
sume a static point source of mass M for the stress-energy 



tensor and solve the field equation for the Schwarzschild- 
like solution. The field equation with this source is 

V^ * = 4:nMS{x) . (2.21) 

Away from the source, the solution to this equation is 



$ = 

r 



(2.22) 



which is, of course, the Newtonian potential. The space- 
time metric is therefore 

The small-$ (small M or large r) limit of the goo compo- 
nent is 500 — > — (1 -I- 2<I>), therefore correspondence with 
Newtonian gravity is assured. 



E. PPN Parameters 

Having found the field equation, one is now compelled 
to ask, how "good" is the theory? We will explore a 
number of definitions of "good" in this paper, but one 
possible definition is to insist that a "good" theory have 
parameterized post-Newtonian (PPN) parameters which 
1) satisfy as many as possible of the solar system gravita- 
tional tests made to date, and 2) match as closely as pos- 
sible the PPN parameters of general relativity. One must 
keep in mind, however, that obviously a theory may be 
"good" in its Newtonian limit yet still give results which 
differ greatly from general relativity in the strong-field 
limit! In fact, we have found no solutions in this SCF 
scalar theory that approximate the cosmological models 
in general relativity. 

Ni H has calculated the PPN parameters for his 
"Lagrangian-based stratified theory", which our simpli- 
fied theory closely resembles. Subtracting the present 
field equation ( |2.7| ) from Ni's shows that the theories dif- 
fer by 



iv^gn.^.'^..-iv^vn,^^'^,. 



-4*\ 



= 4e-^*($o)' + (l-e-^'^)$ 



,00 • (2.24) 

Both of these terms are 0{c^^). Because the PPN formu- 
lation considers terms through 0{c~'^), this theory and 
Ni's must have identical PPN parameters. [The PPN 
approximation is often considered an expansion in pow- 
ers of 1/c since the dimensionless quantities {v /cf' and 
GM /r(? are of comparable size (virial theorem) in bound 
systems, and timelike derivatives do are smaller by a fac- 
tor of order v / c than space derivatives di when the time 
dependence is generated by the motion of sources. The 
Newtonian potential $ ~ —GMjrc? is considered to be 
0(c~^).] The PPN parameters for this theory ^ are 
therefore 



e-7?-0,7-/3 = /3i=/32 = /33-/34 



The theory agrees with general relativity with the ex- 
ception of the Ai parameter. Additionally, due to its 
value of Aijjt disagrees with the Earth rotation rate ex- 
periments. |15[| Essentially, this theory predicts an incor- 
rect relative amount of dragging of inertial frames (goi) 
produced by unit momentum {pv) (see MTW Box 39.2 
[Q). On the basis of PPN parameters, this theory com- 
pares favorably to all of its scalar gravity competitors. In 
fact, one obtains precisely the PPN parameters of Ni's 
SCF theory Q and Rosen's (l6| two field variable theory 
with A = 1, yet both authors postulate a considerably 
more complex field equation than is used here. 



III. THE SCHWARZSCHILD-LIKE SOLUTION 

The PPN formulation predicts how the mid-range field 
of scalar gravity compares to general relativity, but it 
makes no such prediction about the strong-field region. 
A second measure of how "good" the theory is, then, 
is to calculate the analytic results for a simple strong- 
field case and compare those results to the equivalent 
case in general relativity. The metric found in Section 
[ID corresponds to the well-studied Schwarzschild case 



of GR; in this section we will determine the properties of 
the equivalent case for scalar gravity, studying geodesies 
in a manner parallel to §25.3 of MTW |§. 



A. The Effective Potential 



To begin, write (2.22) in spherical coordinates 

ds2 = -e-2M/r-^i2 ^ g2A/A^^2 ^ ^2g2A//r^^2 (3^^) 

The geometry is unaffected by changes in t or 0, therefore 
two constants of motion exist: 



Po == -E , pe = ±L 
Conservation of rest-mass energy, 

gives the orbit equation 



(3.2) 
(3.3) 



- E'^^'l- + e^^'l- (^ ^ + {^\ e-2M/r- + 1 = 

(3.4) 

where A is the path parameter and jf = dr / dX. Defining 
E = -E//i, L = L/ p., A = T / p places the orbit equation 
in its final form 



(2.25) 



%] ^E^-VMrf 



where V^s, the effective potential, is defined as 



(3.5) 




FIG. 1. The effective potential ( p.q ) plotted as a function 
of r /M for a typical specific angular momentum, in this case 
L = 4Af. Notice the "pit in the potential", (see MTW [§, 
Chap. 25]) which gives rise to an innermost stable circular or- 
bit when, for still smaller values of L, the stable and unstable 
circular orbit radii (local min and max of Vdt) coincide, (see 
Section liilBl). 



V,«{rf 



-2M/r 



-2M/r 



(3.6) 



Figure 1 shows a typical plot of the effective potential as 
a function of r. 



B. Innermost Stable Circular Orbit 

One of the first relativistic effects to be noticed was the 
"pit" in the effective potential, the existence of which im- 
plies the presence of an innermost stable circular orbit. 
Any particle that moves within this limit must neces- 
sarily spiral in and impact the singularity at r=0. As 
the present theory also possesses this pit, it is reasonable 
to suppose that it possesses an innermost stable circular 
orbit as well. 

Circular orbits occur at extrema of the effective po- 
tential. Taking the derivative with respect to r of ( |3.6| ) 
returns, after some algebra, the circular orbit angular 
momentum as a function of r: 



^cnc — ' c 



M 



2M 



(3.7) 



It can immediately be seen that no circular orbits exist 
at all below r — 2M . The innermost stable circular orbit 
occurs at the value of r for which Ldrc is a minimum. 
Taking the derivative of (|3.7|) with respect to r yields 



CO = M{i + %/5) 



(3.8) 



This result compares favorably to the predictions of gen- 
eral relativity. GR predicts the innermost stable circular 



orbit occurs at r = 6M, which corresponds to a circum- 
ference of 12'kM. The circumference of the innermost 
stable circular orbit under the scalar gravity theory is 
12.67597rA'f , a difference of only 5.63% 

Well inside the innermost stable circular orbit we can- 
not expect this theory to mimic general relativity. To 
produce the analog of the interior Schwarzschild solu- 
tions would require inclusion of the e~^^ 5ijT^^ term in 
equation 2.7 and the hydrostatic fluid equations, so it 



is possible that the theory would give a maximum neu- 
tron star mass. But since the radial null vectors are 
k± — dt ± e^'^dr, the only possibility for a black hole 
effect would be to find solutions where $ ^ — oo at a 
nonzero value of r. 



IV. GRAVITATIONAL RADIATION 

With the completion of the LIGO gravitational wave 
interferometer nearing, attention in the numerical rel- 
ativity community has been focused upon the produc- 
tion of gravitational wave templates to supply LIGO's 
matched-filtering analysis routines. A third and partic- 
ularly relevant measure of how "good" scalar gravity is 
would be to calculate the frequency of the radiation, as 
observed at infinity, of the innermost stable circular or- 
bit. 

From the scalar gravity Schwarzschild-like metric 

ds^ = -e-^^'/'df + e^^'/'dr^ + r^e^^'/'^ide^ + sin9^dcf>^) 

(4.1) 

the radiation frequency, 2uj, can be defined from the cir- 
cular orbit frequency 



OJ ■ 



d4 ^(^(h ^ _4M/r 

dt dr dt 



L 



pTT ■ (4-2) 



Combining ( p.6\ j and (3/7), yields the effecti ve p otential 
as a function of orbital radius; this value, by (B.5), is also 
equal to the orbital energy squared. The result of this 
algebra is 






- M 



2M 



(4.3) 



A somewhat more useful quantity conceptually, how- 
ever, is the binding energy. 



E^ 



bindint 



i ^circ 



(4.4) 



Substituting ( |4.3| ) and ( p.7[ ) into (O) gives, after some 
algebra, scalar gravity's equivalent of Kepler's Third Law 



Lohl,,e^^''/^-^{r~M)^M 



(4.5) 
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FIG. 2. The frequency spectrum of gravitational radiation 
produced by a test mass in nearly circular orbit about a cen- 
tral object as the test mass adiabatically s piral s inward. The 
graph shows the change in binding energy (4.4) as a function 
of circular orbit frequency (4.2), which is in turn a function of 
the circular orbit radius. The solid line shows the spectrum 
for a test mass moving within a Schwarzschild background 
metric, while the dashed line shows the spectrum for a test 
mass moving within the background metric of ([4.l|). The in- 
nermost stable circular orbit occurs where each curve crosses 
the a;-axis. The agreement between the curves is very good 
almost down to this point. 

Notice that this expression is precisely Kepler's Third 
Law only for large r/M. Although -Ebinding(^) can- 
not be found as an analytic expression, we can eas- 
ily plot the spectrum of the gravitational radiation 
d-Ebinding(w)/da; by knowing both di?binding(w)/dw = 
(di?binding(w)/drcirc)/((iw/(ircirc) and Lu as functions of 
Tcirc- Figure 2 compares the frequency spectrum calcu- 
lated under scalar gravity with that found from general 
relativity. As can be seen, the two theories agree quite 
well almost down to the innermost stable circular orbit. 
The frequency of the innermost stable circular orbit is 
found using ([4.5D evaluated aX r — r^sco ~ ^i{'i + Vs): 



0.0633326 

M 



(4.6) 



General relativity gives a value of Wisco = 0.0680414/Af 
(a difference of 6.92%), again demonstrating the close 
agreement the two theories. This close agreement is not 
too surprising when one considers that there is a sense in 
which the Schwarzschild-like solution of this SCF scalar 
theory "almost" satisfies Einstein's equations. The non- 
zero compo nent s of the Einstein tensor, calculated using 
the metric (4.1) are 



CzQ — -Cz^ — ^2 



G 



-2M/r 



M2 



(4.7) 



rather than zero. This is smaller than the nonzero com- 
ponents of the Riemann tensor in the Schwarzschild met- 
ric by a factor e^^*^/''(Af/r). A dimensionless measure 



of this failure to satisfy Einstein's equations has its max- 
imum at r = M where each nonzero component M^ | G(^ | 
is e^^ = 0.135, from which point they decrease mono- 
tonically. At r = risco this error is 0.0009 and decreases 
rapidly farther out. 



V. DISCUSSION 

The scalar theory presented here is intended as a stand- 
in for general relativity while developing hydrodynnamic 
computer codes for application in problems such as the 
production of gravitational radiation by binary neutron 
star systems. It has the advantage of involving only one 
field variable $ (and its time derivative) a nd s atisfying 
a simple field equation. The field equation (2.7) involves 
only the usual flat spacetime d'Alembertian (nonlineari- 
ties only in terms containing no derivatives) and reduces 
to the flat space linear wave equation outside the region 
occupied by matter. This should simplify its numeri- 
cal treatment, and particularly the imposition of outgo- 
ing wave boundary conditions; there are no curved light 
cones to manage. Its interface to hydrodynamics can be 
identical that of general relativity — it requires a matter 
stress-energy tensor and produces a metric. 

It seems plausible that, in a range including flclds 
stronger than those treated by the PPN methods, this 
theory could also give results that are close to those of 
general relativity. For Coulomb-like gravitational flelds 
our study of the analog of the Schwarzschild metric shows 
results within a few percent of those of GR down to 
the innermost stable circular orbit. For gravitomagnetic 
forces (from metric components g^i generated by mat- 
ter currents pv) it will necessarily fail, as it allows no 
goi metric components. Thus one cannot expect phys- 
ically plausible results when matter velocities approach 
the velocity of light. But since this scalar theory does 
reproduce important qualitative features of general rela- 
tivity (the existence of an innermost stable circular orbit 
for test particles) , it may allow the initial study of impor- 
tant relativistic hydrodynamical phenomena whose accu- 
rate analysis will have to be done using the full Einstein 
theory of gravitation. 
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